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Using a gauge-invariant matrix parametrization of the 
gauge fields, we present an analysis of how the mass gap arises 
in (2+l)-dimensional Yang-Mills theory. We further derive an 
analytical continuum expression for the vacuum wavefunction 
and based on this we calculate the string tension which is in 
excellent agreement with Monte Carlo simulations of the cor- 
responding lattice gauge theory. 



I. INTRODUCTION 

While the perturbative aspects of non-Abelian gauge 
theories were understood many years ago, nonperturba- 
tive phenomena have been difficult to analyze. Detailed 
analyses and calculational techniques are still lacking, 
even though most of the qualitative features are more 
or less clear. 

In particular, the issue of confinement and generation 
of mass gap in Quantum Chromodynamics (QCD) still 
defies a clear quantitative analysis. We have tried to un- 
derstand these nonperturbative mechanisms in a some- 
what simpler but still physically relevant framework, the 
(2+l)-dimensional Yang-Mills theory. To start with we 
have neglected the quark degrees of freedom, keeping only 
the gluons. The confinement mechanism applies to glu- 
ons as well. The lowest excited states of the purely glu- 
onic theory are not the perturbative gluons but rather 
massive, color-singlet bound states of gluons, the so- 
called glueballs. The lowest glueball mass defines essen- 
tially the mass gap of the theory. 

The study of non-Abelian gauge theories in a lower 
dimensional context, such as two spatial dimensions is 
interesting for at least two reasons. Primarily it can be 
a useful guide to the more realistic case of three dimen- 
sions; given the complexity of gauge theories, any further 
progress, even in a simpler lower dimensional context, 
should be interesting. Secondly, gauge theories in two 
spatial dimensions can be an approximation to the high 
temperature phase of QCD with the mass gap serving as 
the magnetic mass. It is the latter connection that makes 
our work relevant in a thermal field theory context, al- 
though there is no explicit mention of temperature. For a 
resume of our work and the relation to (3+l)-dimensioanl 
magnetic mass, see V.P. Nair's contribution to these pro- 
ceedings. 



II. HAMILTONIAN ANALYSIS, 
GAUGE-INVARIANT PARAMETRIZATION 

In an attempt to understand the origin of the mass 
gap and the confinement mechanism, we have analyzed 
SU(N) Yang-Mills theory in two spatial dimensions in a 
Hamiltonian framework in terms of gauge-invariant vari- 
ables [0 One of the advantages of the Hamiltonian 
approach (fixed time) is the fact that the theory is es- 
sentially two-dimensional and various mathematical tech- 
niques from two-dimensional conformal field theory be- 
come applicable and very useful. As is convenient for a 
Hamiltonian formulation, we use the Aq = gauge. The 
gauge potentials are Aj = —it a Af, i = 1,2, where t a are 
hermitian (TV x A)-matrices which form a basis of the Lie 
algebra of SU(N) with [t a , t b ] = if abc + c 
The Hamiltonian can be written as 



c t c , Tr(i a r) = ±<S 
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where e is the coupling constant, Ef = — i-gjtr is the 

electric field and B a = \e ]k {djA k - d k A } + [A j ,A k \) a 
is the magnetic field. In (2+1) dimensions, e 2 has the 
dimension of mass. 

Our basic strategy has been to reformulate the the- 
ory in terms of gauge-invariant variables and analyze the 
corresponding Schrodinger equation directly on the phys- 
ical configuration space C = A/Q (space of gauge poten- 
tials modulo gauge transformations). This is achieved 
by making use of a special parametrization of the gauge 
fields in terms of complex matrices in the following way. 
We combine the spatial coordinates x\, xi into the com- 
plex combinations z — x\ — 1x2, z = x\ + 1x2] corre- 
spondingly we have A = A z = \{Ai + 1^2), A = A Z = 
■^(Ai — iA2) = — {A z y . The parametrization we use is 
given by 



A z = -d z MM~ 



A- z = M t " 1 9jM t (2) 



Here M, are complex SL(N, C)-matrices. 

From (||) it is clear that the gauge transformation, 
A. L — > gAig^ 1 — digg" 1 , is expressed in terms of M, 
by 



M -> gM , 



.9 



(3) 



for g(x) £ SU(N). In particular, if we split M into a 
unitary part U and a hermitian part p as M — Up, then 



1 



U is the 'gauge part' and it can be removed by a gauge 
transformation. The combination H — M'M = p 2 is 
a gauge-invariant field and it provides a parametrization 
for the configuration space C. The hermitian field H 
belongs to SL(N,C)/SU(N). 

The parametrization (^) introduces a new symmetry in 
the theory, what we call "holomorphic" invariance. This 
arises as follows. The original gauge fields A, A remain 
invariant under the transformation M — > MV(z), M' — > 
V{z)M\ H ->■ V{z)HV(z), where V, V depend only on 
the holomorphic coordinate z and the antiholomorphic 
coordinate z respectively. Since the original theory is 
unaffected by this transformation, we should require that 
all physical states remain invariant under "holomorphic" 
tranformations . 

Having derived a particular parametrization for C one 
is able to ask more detailed questions regarding its struc- 
ture. An important question is: what is the integration 
measure (volume element) on C ? One of the advantages 
of the particular parametrization that we introduced is 
that we can give an explicit answer to the above question. 

The volume element in the original space of gauge po- 
tentials A is given by 

dn{A) = Y[dA a (x)dA a {x) = \det(DD)\ dfi(M,M^) (4) 



the physical two-dimensional space corresponding to the 
coordinates z, z. The integrand thus requires an exten- 
sion of the matrix field H into the interior of M 3 , but 
physical results do not depend on how this extension is 
done H . Actually for the special case of hermitian ma- 
trices, the second term can also be written as an integral 
over spatial coordinates only ||. 

The inner product for gauge-invariant physical states 
is given by 



(1|2> 



d^(H)e 2cAS{H) V1(H)V 2 (H) 



(9) 



Here we begin to see how conformal field theory could 
be useful; equation (||) shows that matrix elements of 
the (2+l)-dimensional gauge theory can be thought of as 
correlators of the hermitian WZW model, which is a con- 
formal field theory. This point of view allows us to derive 
further constraints on the dependence of the wavefunc- 
tion on the gauge-invariant variables H. The correlators 
of the hermitian WZW model have been studied before 
H and following conformal field theoretic arguments one 
can show that all wavefunctions of finite inner product 
are functions of the current 



J a = —(dHH~ 1 ) a 



(10) 



where D, D are covariant derivatives in the adjoint rep- 
resentation; D = d - dMM-\ D = 5 + Aft-i<9Mt. 
dp(M, M') is the standard volume element in the space 
of SL(N, C)-valued fields. Factorising M into a unitary 
part U and a hermitian part p, M = Up, dp(M, M<) can 
be written as 



dp(M, M f ) = dp(H)dp(U) 



(5) 



where dp(U) is the Haar measure for SU(N) and dp(H) 
the corresponding measure for hermitian fields. The vol- 
ume element for C is now obtained as 

MC) = J^) = ^ et{DD)mH) (6) 

The problem is thus reduced to calculating the deter- 
minant of the two-dimensional operator DD. Using a 
gauge-invariant regulator we find 0,p[ 



dn(C) - e 2cAS{H) dp{H) 



(7) 



where ca is the adjoint quadratic Casimir (ca = N for 
SU(N)) and S{H) is the Wess-Zumino-Witten (WZW) 
action for the hermitian matrix field H given by 

S(H) = — Tr / d 2 x dHdH' 1 
2tt J 

+ ^-T r y d 3 xe^ a H- 1 d tl HH- 1 d u HH- 1 d a H (8) 

As is typical for the WZW action, the second integral 
is over a three-dimensional space M 3 whose boundary is 



The result is not surprising. The Wilson loop operator 
itself can be written in terms of the current as 



W{C) 



Tr P e 



— <t (Adz+Adz) 



Tr P e 



Since in principle all gauge-invariant functions of A, A 
can be constructed from W(C), the current J should 
suffice to generate all the gauge-invariant states of the 
theory. It is thus important to express the original Hamil- 
tonian in terms of the current J, which can be thought 
of as a collective field. 

The potential energy term can be written as 



where m — !L ^ A . 



mcA 



dJ a (x)dJ a {x) (12) 



The kinetic energy term becomes 



r{x) 



VL a b{x,y) 



5J a (x) 

CA S ab 

7r 2 (x — y) 2 



6J a (x) 5J b {y) 



. fab C J c (y) 
tt(x - y) 



(13) 



The first term in T essentially counts the number of J's 
and the second term replaces pairs of J's by the operator 
product expansion for currents in the WZW-model. It is 
the first term which is responsible for the appearance of 
a mass gap. As it is obvious from (|l^), the lowest excited 
eigenstate of T is the current J itself: 



TJ a (x) = mJ a (x) 



(14) 
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The current J can be thought of as the nonperturba- 
tive gauge-invariant gluon with dynamical mass to. Al- 
though J a is an eigenstate of T, it is not an accept- 
able physical state. As we mentioned earlier, physical 
states should be invariant under "holomorphic" transfor- 
mations. The state with the minimum number of J's 
which satisfies "holomorphic" invariance is a 2 J-state of 
the form dJ a dJ a . This is also a color-singlet state as 
expected. Such an expression after a careful regulariza- 
tion ("holomorphic" invariant point splitting) is a good 
candidate for the lowest glueball state ++ 

$2 = J f(x,y)Tt : dJ(x)(H(x,y)H~\y,y))dJ(y) : (15) 

Such a state has a mass gap of at least 2m. The regulator 
term introduces interactions among currents which make 
the estimation of the binding energy rather difficult. 



III. VACUUM WAVEFUNCTION AND STRING 
TENSION 

So far we have considered the spectrum of the kinetic 
energy term T. Notice that the potential energy term 
V, (p^), gives contributions of the order k 2 /e 2 ~ k 2 /m, 
where k is a typical momentum. For momentum modes 
k -C e 2 ~ to, the kinetic energy term is dominant and its 
spectrum provides an approximation for the spectrum of 
the whole theory. In this limit, which can be thought of 
as the strong coupling limit, the potential energy term V 
can be treated in a perturbative fashion. Let us see how 
this works in the calculation of the vacuum wavefunction. 

As far as the kinetic energy operator is concerned, 
$0 = 1 may be taken as the vacuum wavefunction. 
This might seem rather trivial, however it is important 
to notice that the corresponding norm ($o|^o)j which 
is also the volume of the configuration space C, coin- 
cides with the partition function of the hermitian WZW- 
model, which upon appropriate regularization is finite Q . 
In other words $o is normalizable with the inner product 
(0) . Essentially the exponential factor in (|9|) provides the 
necessary damping to render the norm finite. (It is infi- 
nite in the Abelian case ca = 0, where the exponential 
factor is absent.) 

For low momentum modes, k <C to , the inclusion of 
the potential energy term leads to a modified vacuum 
wavefunction which can be written as 



(16) 



where P is a functional of the J's which can be expanded 
in powers of 1/to. The various terms in this expansion 
can be determined from the Schrodinger equation 



Wo = (t + v)v q = 

The first few terms are given as H 



(17) 



P = - 



m 2 CA 



Tr / : dJdJ 



^-) 2 Tr J [: dJ{Vd)dJ+ l -dJ[J,d 2 J] :] 

- 2 ( Tr [ \: dJCDdfdJ +l\VdJ, dJ]d 2 J 
\m 2 c A J J 9 

+ ^[v5 2 j, J]d 2 j-\\J, Bj][Bj, B 2 j] 

9 6 
-^[J,dJUd 3 J] :] + 0(±) (18) 

where T>h = ^dh — [J, h]. The normal ordering of vari- 
ous terms in ( p^ ) is necessary for P to satisfy (|l7|). The 
second derivative term in ( |l3| ) can give singularities when 
acting on composite operators. The normal ordering sub- 
tracts out precisely these singularities. 

Our approach provides the first analytical, continuum 
calculation of the vacuum wavefunction. There are sev- 
eral interesting points regarding this. The leading order 
term for the vacuum wavefunction is 



"to w exp 



exp 



m 2 CA 



Tr / : dJdJ 



2me 



rTr / B 



(19) 



The calculation of expectation values involves averag- 



ing with the factor \&o*o 



where S, as seen 



from the above equation, is the action of a Euclidean 
two-dimensional Yang-Mills theory of coupling constant 
g 2 = me 2 — e^CA/^ir. Thus, retaining only the leading 
term in 4"0) the expectation value of the Wilson loop in 
the fundamental representation is given by 



(W F (C)) = exp 



e cacf 

4-7T 



(20) 



where Ac is the area of the loop C and cp is the quadratic 
Casimir of the fundamental representation. The expec- 
tation value of the Wilson loop exhibits an area law be- 
havior, as expected for a confining theory. From ( p0| ) we 
can identify the expression for the string tension a as 



e 4 c A c F ±fN 2 -l 

— e 



4tt 



8tt 



(21) 



Recent Monte Carlo calculations of the string tension 
give the values V^/e 2 = 0.335, 0.553, 0.758, 0.966 for 
the gauge groups SU{2), SU(3), SU(4) and SU{5) re- 
spectively H . The corresponding values calculated from 
d|) are 0.345, 0.564, 0.772, 0.977. We see that there 
is excellent agreement (upto ~ 3%) between (|2l]) and 
the Monte Carlo results. It is further interesting to no- 
tice that our analytic expression for the string tension 
( pl| ) has the appropriate A^-dependence as expected from 
large- A calculations. 

Another interesting point is the fact that we have ob- 
tained an expression (O) local in terms of the current J 
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but nonlocal in terms of the original magnetic field B. 
All previous attempts to derive an analytical expression 
of the vacuum wavefunction are based on the assumption 
that, in the strong coupling regime, such a wavefunction 
can be written as a local expansion in B |^,[n)]]. Our 
results obviously do not agree with this approach. It is 
thus interesting to investigate whether one could incorpo- 
rate these nonlocal terms and verify our results in lattice 
calculations. 

The approximation of the vacuum wavefunction by the 
first few terms in (|l8|) makes sense only in the low mo- 
mentum regime k -C m. Although very cumbersome, in 
principle we could determine arbitrarily many terms in 
the series expansion ([l8]), so one would expect that by 
resumming up various terms one might be able to get 
information on the vacuum wavefunction away from the 
low momentum region. In particular, what about the 
short distance behavior, k ^> e 2 ? Indeed this can be 
retrieved by converting the l/m expansion into a series 
expansion in J's. The terms in (|l8|) can be naturally re- 
arranged into terms with two J's, terms with three J's, 
etc. The series of terms (infinite many) with only two 
J's can be summed up to give 



P = 



1 



B a {x) 



x.y 



1 



(to + \J m 2 — V 2 ) 



B a (y) (22) 



(T + V)J a e p = yjm 2 - V 2 J a e p 



(24) 



where in the expression for P we have kept only the lead- 
ing 2J-term contribution and have neglected the higher 
J-terms. As expected in a relativistic theory the mass 



gets corrected to its relativistic expression yf 



Jc 2 



which also displays the correct high and low momentum 
limits. In a similar fashion we can obtain expressions in- 
terpolating between short and long distance regimes for 
other quantities of interest. 
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In the weak coupling or high momentum regime J:>e 2 , 
the leading order term in ( |22| ) is 



p = ~I nj,) 



x.y 



(-V 2 



B a (y) (23) 



x,y 



This is the vacuum wavefunction for an Abelian theory 
as expected. Of course, in the low momentum regime 
k <C e 2 the leading order term in ( |22] ) will reproduce 

Now what about the 3 J-terms and higher order terms? 
Part of the 3 J-contribution is just what is needed to co- 
variantize the Laplacian operator V 2 in (p2|), but that is 
not all. We can explicitly derive the full 3 J-contribution 
What is important is the fact that comparing the 
2 J and 3 J-terms we find [|J that the 3 J-term is subdom- 
inant compared to the leading 2J-term in P for both 
the low and high momentum regimes. Similar arguments 
hold, based on dimensional analysis, for the higher J- 
terms. So the expansion of P as a series expansion in J's 
is consistent with both low and high momentum regimes. 

So far we have discussed the structure of the vacuum 
wavefunction. One could, in principle, extend the above 
analysis for the excitation spectrum of the Schrodinger 
equation. For example, in the absence of the potential 
term V, the current J is an eigenstate of the kinetic en- 
ergy operator T, (|l4|), with eigenvalue to. One could 
now ask what the corresponding modified eigenstate and 
eigenvalue should be, once the potential term is included. 
We find 
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